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Abstract 

We use a boost-invariant one-dimensional (cylindrically symmetric) fluid dy- 
namics code to calculate thermal photon production in the central rapidity 
region of S+Au and Pb+Pb collisions at SPS energy (i/s = 20 GeV/nucleon). 
We assume that the hot matter is in thermal equilibrium throughout the ex- 
pansion, but consider deviations from chemical equilibrium in the high tem- 
perature (deconfined) phase. We use equations of state with a first-order 
phase transition between a massless pion gas and quark gluon plasma, with 
transition temperatures in the range 150 < T c < 200 MeV. 
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I. INTRODUCTION 



Electromagnetic probes (photons and leptons) have long mean free paths in hadronic 
matter and are therefore well suited for studying the earliest stages of ultrarelativistic nuclear 
collisions. Preliminary single photon spectra from 200 GeV/nucleon S+Au collisions at 
CERN's Super Proton Synchrotron (SPS) have recently become available from the WA80 
experiment, Q and results from Pb+Pb collisions should be available in the near future. 
In the present work we calculate thermal photon spectra in the central rapidity region for 
SPS collisions, based on a boost-invariant fluid-dynamical description of the hot matter; a 
similar approach was used in Ref. [Q]. 

Earlier predictions of photon spectra f3[ and of dileptons from several resonances Q || 
show a useful sensitivity to the assumed QCD transition temperature when the equation of 
state incorporates a strong first-order phase transition and transverse expansion is neglected. 
The correlation between the apparent temperatures of photon spectra and the assumed 
transition temperature was suggested as a thermometer for the transition temperature. || 
One question we wish to study here is whether a useful correlation can still be found if 
transverse expansion is taken into account. 

In our calculation we assume a boost-invariant longitudinal expansion as discussed by 
Bjorken, |7J coupled to a cylindrically symmetric transverse expansion. The one-dimensional 
fluid-dynamics code we developed uses the Godunov method as described by Blaizot and 
Ollitrault. || We assume thermal equilibrium throughout the evolution of the system, but 
consider deviations from chemical equilibrium in the high-temperature phase by allowing 
the quark and antiquark densities to be a (fixed) fraction of the equilibrium value. In 
the photon production rate we include the effect of the a>i resonance. We calculate the 
transverse momentum distribution for photons in the range 1 < px < 2 GeV. We investigate 
the sensitivity to different assumptions about the initial temperature, freezeout temperature, 
and quark fraction, and compare production rates to the preliminary WA80 data. We use 
standard high-energy conventions, c = h = ks = 1- 



II. FLUID-DYNAMICAL EVOLUTION 

For high collision energy we expect longitudinal boost invariance, so the behavior 
of the produced matter at different rapidities is the same in the longitudinally comoving 
frame for fixed proper time r = \ft 2 — z 2 , where z is the distance along the beam axis. At 
t = the colliding nuclei reach the point of maximum overlap and are assumed to form a 
longitudinally expanding pancake. The hot matter has thermalized at r = tq (^ 0.2 fm/c 
|p|JT0|1 ) , when transverse expansion begins, coupled to the longitudinal expansion. We treat 
the hot matter as one- dimensional, assuming cylindrical symmetry and boost invariance 
with equations of motion || 

A(rrT 00 ) + |-(rrT r0 ) = -rP, (la) 
or or 

^( rrT ° r ) + |:( rrTrT )= rR ( lb ) 
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Here P is the pressure, r is the radial coordinate, and the energy-momentum tensor T^ v is 
jifxv = ^ e _|_ p^n^ii u — P g^ u , with e being the energy density, it the four- velocity, and metric 
tensor g^ v = diag(l, — 1, — 1, — 1) . The r— weighting takes into account the radial geome- 
try, and the r— dependence is due to the coupling between the longitudinal and transverse 
expansions. 

The fluid dynamics code uses the Godunov method as described by Blaizot and Olli- 
trault. || The system is described by the relativistically covariant energy-momentum tensor 
elements T 00 and T 0r . The program calculates the changes in the cell averages of these 
tensor elements, respectively 00 and O 0r , by calculating the flows T 0r and T rr across the 
cell walls. A relationship between the flows at the walls and the changes in Q^ v is found 
by integrating (Jl|) over each space-time cell, assuming plane similarity forms for the flow 
patterns. We approximate the space-time averaged pressure (P), by using the initial value 
of T^ v to estimate at intermediate times, and constructing (P) from these estimates for 
00 and Or . This estimated value of (P) is then used to solve for Q^ u at the end of the 
time step. At the beginning of the next time step, we set T 00 and T 0r throughout each cell 
equal to O 00 and O 0r for that cell at the end of the previous step. 

From r = until the transverse expansion starts at r = r , we assume a boost-invariant 
cylinder of radius R < (the radius of the smaller nucleus), filled uniformly with QGP at 
temperature T = T . This is approximately compatible with the initial entropy density for 



short times. [11 1 We determine Tq by assuming entropy conservation for r > tq, hence 



s(T 



3.6 dN/dy 
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where s is the entropy density, with total (charged plus neutral) multiplicity density dN/ dy 
240 and 1600 for S+Au and Pb+Pb collisions, respectively. 

The equations of state (EOS's) that we use here are of the form 
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where g q is the number of massless degrees of freedom in the deconfined phase. We treat 
only the case of zero baryon density, so the entropy density is 
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(4) 



T 



independent of the phase of the matter. Below T c , the EOS is that of a massless pion gas. 
Because recent calculations have predicted that the quarks may reach only a fraction of their 
equilibrium number by the beginning of transverse expansion we take g q = 16 + 21s, 



where x is a parameter that we vary between and 1 to simulate the effect of reducing the 
quark density in the QGP below the equilibrium value (x=l is equilibrium for two flavors 
of massless quarks). The vacuum energy density in the deconfmed phase, B, is related to 
the transition temperature, T c , by requiring equal pressures in the deconfmed and hadronic 
phases at T = T c : 

B = (5) 
90 V ; 

For all the fluid dynamic calculations we have taken the radial cell size 5r and time 
step 5t small enough to be within a few percent of the continuum limit. We also required 
5r <C Sr/2, guaranteeing that a similarity pattern from one side of a cell does not overlap 
the opposite side of the cell or the pattern from the other side. Also, because the radial 
pattern size is much smaller than 8r, we may assume plane similarity solutions as a good 
approximation of the intercell flow in spite of the radial geometry of the system. 



III. THERMAL PHOTONS 



The central region photon p? distribution from boost-invariant hot matter is 
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where 
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sin 2 9 + 2v sin 9 cos (J) + v 2 (l — sin 2 9 sin 2 



(6) 



(7) 



Here R is the photon production rate per unit four-volume in the CM frame of the hot 
matter, and v is the transverse velocity of the matter in the cell characterized by proper 
time r, space-time rapidity rj and radial position r [measured in the frame moving with 
transverse velocity zero and longitudinal velocity tanh rj in the lab]. Evaluating the integrals 
over rj and p with the (^-functions, we obtain 



d 2 AL 



p T dp T dy 



p* 



drr I dr2irr / d9 sin 9 / u Y > i ,., 

o Jo p^ \d 6 p / 



p=p* 



Pt 

7 sin 2 9 + 2v sin 9 cos + v 2 ( 1 — sin 2 9 sin 



1/2' 



(8) 
(9) 



4 



The photon production rate from thermally and chemically equilibrated quark-gluon 
plasma is 0] 



where a = 1/137, and a s = c/ 2 / 47r = 0.4. [We use the semi-empirical formula of Ref. ||12| , 
which is almost identical to the exact rate.] 

Equation fllOf) is first order in a s (or g 2 ), and as g is of order unity, its accuracy is of the 
order of unity. We also use it for the hadron gas, as its uncertainty is larger than the difference 
between the first-order QGP and hadron gas production rates. ]13| The contribution of the 



a\ resonance as an intermediate state was discussed recently by several authors. [13,15,^] 
Here we include the contribution of the a± meson as parametrized in Ref. ||14|| . 

The calculation follows the space-time evolution of the collision based on eq. (|8|) and 
results in nearly-exponential spectra. We then mimic the experimental procedure and extract 
a temperature (the fit temperature, Tf it ) from the obtained spectra. For fitting purposes, 
we take (|8]) with v = (i.e. ignoring transverse expansion), resulting in the formula || 



p T dp T dy ^ 



which is accurate to about one percent for the range we consider here. 

The photon production in QGP is dominated by quark-gluon collisions. To simulate 
scenarios where the quark number is smaller than the equilibrium value, we multiply the 
rate (|10D by x, the fraction of the equilibrium number of quarks, for the QGP component of 
the fluid. This representation of the deviations from chemical equilibrium is valid whenever 
the quarks and antiquarks are non-degenerate, so that Maxwell-Boltzmann statistics can be 
used in place of Fermi statistics, and when the production from quark- ant iquark annihilation 
(which is proportional to x 2 ) is small. 



IV. RESULTS 

Our standard calculation uses tq= 0.2 fm/c, x=l, freezeout temperature Tf D = 100 MeV, 
and includes the a± resonance; if a different value of one of the parameters is given, it is to 
be assumed that the others are held at the standard values. In Fig. 1, we vary T c from 150 
to 200 MeV and calculate the resulting values of the temperature, Tf it , needed in eq. ( |TTD to 
fit the resulting spectrum for a central S+Au collision at SPS energy. Our standard value 
for the equilibration time, tq= 0.2 fm/c, implies an initial temperature of To = 355 MeV for 
the S+Au system. 

In the absence of transverse expansion, Tf it is a monotonically increasing function of T c , 
so one can infer T c given Tf it from the measured photon spectrum. [0] This can be seen in 
the lowest curve, where we have run our simulation without transverse expansion. 

Including transverse expansion, we find that we can no longer infer T c from Tfu, due to 
the non-monotonicity of the relation (except for long thermalization times; see the curve 
with r =1.0 fm/c, corresponding to T = 208 MeV). This demonstrates the large effect of 
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the transverse motion of the fluid on the photon spectrum: lowering T c decreases the bag 
constant, which increases the ratio P/e in the QGP, producing more transverse motion. The 
increase in caused by the fluid motion is enough to overcome the decrease of Tf# with 
decreasing T c found without transverse expansion. This non-monotonicity occurred for all 
parameter sets we tried, although for some (e.g., tq = 1 fm/c) the minimum appears outside 
the range of T c investigated here. As the initial temperature is lowered, the value of T c below 
which the effects of flow dominate also decreases. This may be attributed to the fact that, 
as T approaches T c , less flow develops before the mixed phase is reached. 

Contrary to experience without transverse expansion, ]3[ increasing Tf a to 150 MeV 
decreases Tf it , another indication that the transverse motion has a large effect on the spec- 
trum. Apparently the lower-temperature hadronic matter expands outward at a velocity 
high enough that the contribution to px from this motion more than makes up for the low 
temperature. Thus, removing the contribution from the late stages (by raising Tf ) decreases 



For i / 1, the quantity we call T c on the figure is the chemical equilibrium (x — 1) 
transition temperature, 



rather than the larger non-equilibrium transition temperature used in the simulation, which 
is given by eq. fl5|). We do this because the equilibrium transition temperature is more useful 
for comparisons with theoretical calculations. The curve for x = 0.9 shows that Tju is about 
3 MeV higher than for x — 1. The increase occurs because the temperature of the mixed 
phase, which makes a large contribution to the photon spectrum, is raised when x is lowered. 
We do not show results for other values of x; the change in Tfa is approximately linear in x. 

The effect of increasing tq to 1 fm/c is dramatic: Tfu becomes monotonic within the 
range of interest, and is lower than for the other parameter sets considered. As mentioned 
above, this is due to the decrease in the role of the transverse expansion, which is assumed 
to start at r . 

Figure 2 displays the effect of varying the parameters on the spectra, keeping the value 
of the critical temperature fixed (T c = 170 MeV). We plot d 2 iV 7 / pTdprdy vs. px for different 
scenarios. The spectra are nearly exponential. It is seen that increasing the equilibration 
time (i.e. decreasing the initial temperature) decreases the high p? contribution to the pho- 
ton spectrum and thus decreases the fit temperature. Increasing the freezeout temperature 
removes a fraction of the radiated photons and decreases the yield, also as expected. Ex- 
cluding the a\ resonance decreases the yield by about one-third in the transverse-momentum 
range investigated. This is roughly consistent with several other calculations. [|T^-[T6|j2]] 

Figure 3 shows d 2 N 1 /pTdpxdy vs. px (per central collision) for the preliminary WA80 
data (assuming a minimum bias cross-section of 3600 mb), along with calculated spectra 
using selected parameter sets. The data between 1 and 2 GeV are best fitted with Tf it = 213 
MeV, using eq. (ITT]) . We use this value of Tf it and Figure 1 to extract the value of T c (between 
150 and 200 MeV) for each parameter set. The experimentally determined fit temperature, 
Tf it — 213 MeV is below the values predicted by the model for some parameter choices. In 
these cases T c was estimated by minimizing Tf it . 



fit- 




(12) 
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For the standard parameter set and for x = 0.9 we used T c =170 MeV. Note, however, that 
for i^l, the non-equilibrium transition temperature is larger than this, as explained earlier. 
The difference is sufficient to counteract the decrease in the number of degrees of freedom 
in the plasma, leading to a somewhat higher yield than with the standard parameters. Due 
to the large error bars of the preliminary WA80 data, the calculated spectra are consistent 
with the data points, although generally below them for these parameters. In the absence 
of transverse expansion, photon production is increased by approximately a factor of two by 
including pion masses, || so agreement with the preliminary WA80 data will probably be 
further improved with a more realistic equation of state. 

Increasing the equilibration time to Tq—1 fm/c (and accordingly decreasing the initial 
temperature) leads to T c =190 MeV. The increase in T c increases the photon yield. In 
addition, the lower initial temperature and the higher critical temperature results in less 
flow being built up by the time the system reaches the mixed phase than with the standard 
parameter set. Therefore, the system spends a longer time in the mixed phase for this set 



of parameters. This also increases the yield, as pointed out by Shuryak and Xiong. [16 
To quantify the difference between these scenarios, we "measured" (in the code) the total 
four- volume associated with the mixed phase and found it to be about 25% larger for tq—1 
fm/c than for r =0.2 fm/c . Our results with the long equilibration time (t =1 fm/c) tend 
to overestimate the data. Excess photons come from both the QGP and the hadronic phases 
in our calculation relative to Ref. Jl6| . 

Increasing the freezeout temperature to T/ o =150 MeV decreases the photon yield partly 
because a contribution is removed, and partly because we use T c = 160 MeV (rather than 
170 MeV) in this case. Note that 160 MeV< T c < 170 MeV all give approximately the 
same Tf it for Ty o =150 MeV, though the photon spectrum has a larger magnitude when T c 
is increased due to the large contribution from the mixed phase. Thus, even if T c cannot 
be determined by Tju alone, it may be possible to ascertain T c by comparing to the total 
production rate as well. 

Figure 4 displays predicted spectra for a central Pb+Pb collision at SPS energy, using 
the parameter sets from Fig. 3. The initial temperatures corresponding to r = 0.2 fm/c and 
t = 1 fm/c are T = 441 MeV and T =258 MeV, respectively. The shapes of the spectra and 
the relative production rates are similar to those for the S+Au collisions. For comparison, 
we also show the spectrum obtained in the absence of transverse expansion, taking T c = 200 
MeV. We find that transverse expansion has a large effect on the production rate, in apparent 
disagreement with Ref. j2|. 

Our conclusion is that transverse expansion destroys the correlation suggested in Ref. [|3| , 
so that T c cannot be determined just from measurements of the shape of the single photon 
spectrum, as parametrized by Tf it . The monotonic dependence of Tf it on T c is destroyed 
because the transverse flow increases as T c decreases. However, the total production rate 
is also sensitive to T c , so it is possible that comparison to the production rate instead of 
just the shape of the transverse momentum spectrum may be sufficient to determine T c from 
data. In the transverse-momentum range studied, the inclusion of the ai resonance increases 
the photon yields by about 50%. Preliminary WA80 data appear to rule out large values 
for the equilibration time. 
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FIGURE CAPTIONS 



Fig. 1: Tfn vs. T c for central S+Au collisions at SPS energy with various parameter sets. 
The standard set consists of r = 0.2 fm/c, (corresponding to T = 355 MeV), x—1, 
Tj = 100 MeV, and includes the a± resonance; the other sets differ from the standard 
one in the value of the parameter indicated. 

Fig. 2: Photon spectra with different parameter sets keeping T c = 170 MeV constant. The 
standard set is the same as in Fig. 1. 

Fig. 3: The preliminary WA80 data |1| and selected fits. For each curve the value of T c 
giving Tf it =213 MeV from Fig. 1 is used as explained in the text. 

Fig. 4: Predicted photon spectra for central Pb+Pb collisions at SPS energy. 
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This figure "figl-l.png" is available in "png" format from: 
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This figure "figl-2.png" is available in "png" format from: 
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figure 3 
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This figure "figl-3.png" is available in "png" format from: 
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figure 4 
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